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DIVEEGENCE OF SWEPT WINGS 
By Franklin W. Diederich and Bernard Budiansky 


SUMMARY 


An analyais of the divergence of swept untapered and tapered wings 
with stiffnesses varying as the fourth power of the chord has been 
performed and checked experimentally. The results are presented in a 
set of charts and approximate formulas suitable for quick estimates of 
the divergence dynamic pressure and hence the divergence speed. 

.These results indicate that the divergence speed drops rapidly 
as sweepforward increases to about t0° but that wings with moderate or 
large sweepback cannot diverge. The location of the elastic axis is 
found to affect the divergence speed most at low angles of sweep, where 
movement of the elastic axis forward (or the aerodynamic center aft) 
raises the divergence speed. The effect of wing taper is to increase 
the divergence speed of essentially unswept wings and to decrease the 
divergence speed of wings with moderate and large angles of sweep in 
the case of the prescribed stiffness variation. Evidence is presented 
to indicate that these effects may not be observed for act ua l stiffness 
variations, in which cases a more refined analysis must be resorted to. 


INTRODUCTION 


The emphasis on the use of sweptback or sweptforward wings for 
high-speed flight has created widespread interest in the aeroelastic 
behavior of swept wings. The present paper is concerned with the 
theoretical determination of one of the most fundamental aeroelastic 
parameters, the wing divergence speed.. 

The' divergence of wings or tail surfaces is an instability 
phenomenon which results from the interaction of aerodynamic and 
structural forces. If a wing or tail is given a deflection of 
arbitrary magnitude, the aerodynamic forces often act in such a way 
as to increase the given deflection, whereas the structural forces 
always tend to decrease the deflection. Since the aerodynamic forces 
increase with the flying speed, whereas the structural forces are 
independent of it, a speed will often exist at which the two sets of 
forces are exactly in balance, so that they tend to maintain the given 
deflection. This speed is known as the divergence speed, since any 
further increase in speed causes the aerodynamic forces to predominate 
over the restraining structural forces and tends to increase any 
deformation until structural failure occurs. 




4 


2 



The theory of divergence of unswept wings has reached a considerable 
degree of refinement. Since H. Reissner's original analysis of the 
problem (reference l), a great deal of material has been published on 
this subject both in the United States and in Great Britain (for instance, 
references 2, 3, and 4); the latest methods treat arbitrary stiffness 
variation (reference 2) and account for aerodynamic induction 
(references 3 and 4). 


The analysis of the divergence of swept wings is complicated by 
the fact that, unlike the case of unswept wings, the air forces depend 
on the bending deformations as well as on the twisting deformations. 

A first approximation to the solution of the problem is presented in 
reference 5 by means of a "semirigid" approach that does not take into 
account spanwise variation of the wing distortion. 

The present paper gives the results of a more exact analysis that 
considers the effect of elastic bending and twisting along the wing span. 
Theoretical derivations are given for the divergence speeds of swept 
wings with constant chord and constant flexural and torsional stiffnesses 
along the span, as well as for swept wings with linearly varying chords 
and with flexural and torsional stiffnesses varying quartically with 
the chord. The results of the theoretical analysis are presented in 
curves of nondimensional parameters from which the divergence speed can 
be estimated for a given design. 


In order to verify the theory and establish the effects of the 
assumptions involved, a limited number of tests were made in the 
Langley 4.5—foot flutter research tunnel'on models of constant chord 
and stiffness at low Mach numbers. The results of these tests are 
presented and are compared with the theory. 



SYMBOLS 


aspect ratio (=■ S? 311 - 

VTotal area 


dimensionless parameter 


q cos^\m e e^c r 2 L 2 
(GJ)„ ~~ 


value of parameter a at divergence 

dimensionless parameter 

value of parameter d at divergence 

chord, measured perpendicular to elastic axis, feet 


q cos 2 Am e c r L3tan A 
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chord at effective root (fig. l), feet 
effective chord at tip (fig. l), feet 

effective tending stiffness in planes perpendicular to 
elastic axis, pound—feet^ 

bending stiffness at effective root, pound—feet^ 

distance from elastic axis to aerodynamic center 
(positive forward), fraction of chord 

effective torsional stiffness in planes perpendicular to 
elastic axis, pound-feet^ 

torsional stiffness at effective root, pound—feet^ 

chord ratio, c/c r 

constants 

running air load along elastic axis, pounds per foot 
length of one wing along elastic axis (see fig. l), feet 
free—stream Mach number 

critical Mach number of section perpendicular to 
elastic axis 

section lift—curve slope, per radian 


effective section lift—curve slope, per radian 

dynamic pressure, pounds per square foot 

dynamic pressure at divergence, pounc 3 per square foot 


At d 2 

2 D . 


free—stream velocity, feet per second 

free—stream velocity at divergence, feet per second 

component of free-stream velocity normal to elastic axis, 
feet per second 

distance along elastic axis (see fig. l), feet 



a. 


T) 

r 

A 

\ 

p 


effective angle of attack of a section perpendicular to 
elastic axis, radians 

semispan position along elastic axis (j/L) 

local dihedral angle, radians, or slope of wing deflection 
curve at elastic axis (see fig. l) 


angle of sweep measured to elastic axis, positive for 
sweepback, degrees 


taper ratio 



density of air stream, slugs per cubic foot 
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angle of twist in a plane perpendicular to elastic axis, 
radians ' ’ 


THEORETICAL RESULTS 
Untapered Swept Wings 


A theoretical analysis of the divergence of a swept wing of uniform 
chord (fig, 1(a)) and stiffness is contained in the appendix. The 
analysis involves the following limitations and assumptions: 

(a) Aerodynamic induction is taken into account only insofar as 
an over-all correction is applied to strip theory 

(b) Aerodynamic as well as elastic forces are based on the 
assumption of small deflections 

(c) The wing is clamped at the root perpendicular to a straight 
elastic axis (see fig. 1 (a)), and all deformations are given by the 
elementary theories of bending and of torsion about the elastic axis. 

Within the limitations of these assumptions an exact theoretical 
solution for the divergence speed is obtained. The solution consists 
of a relationship between two nondimensional parameters, 

q. D cos 2 Am e e 2 C r 2 L 2 

D ( n 


( 1 ) 


( 2 ) 


d D = 


cos 2 Ajn e c r L3tan A 

(El)r 


which ia presented in figure 2. 

It i3 seen that the ratio of these parameters 

a D 


(GJ) r 

f 

tr< 

_(El) r _ 

e l c rJ 


tan A 


(3) 


depends only on the known geometrical and physical parameters of -the 
problem. Thus, from the theoretical plot of against dp/ap shown 

in figure 3> 'the divergence speed can he obtained for any particular 
uniform wing. In order to cover the entire range of values of the 
independent parameter, it is convenient to .plot dp against ap/dp for 

large positive and negative values of dp/ap. 


Tapered Swept Wing 3 

A theoretical solution for the divergence speed of linearly tapered 
swept wings (see fig. 1(b)) has also been effected in the appendix. 

The assumptions listed in the foregoing section apply to this case as 
well. The bending and torsional stiffnesses are assumed to vary as 
the fourth power of the chord; this variation is realized for a wing 
having geometrically similar cross sections, such as a solid wing or one 
all structural dimensions of which are proportional to the local chord. 

Theoretical curves similar to those of figure 3, relating the 
nondimensional parameters a D and dp, have been computed for swept 

wings having taper ratios \ of 0.2, 0.5, and 1.5. The various branches 
of these curves, together with those for the uniform wing (taper 
ratio 1.0) are, for convenience, given separately in figure 4(a), 4(b), 
and 4(c). 


Effective Lift-Curve Slope 

In the calculations of the lift acting on the individual wing sections 
it has been assumed that the aerodynamic interaction of the sections may be 
neglected if an over-all span correction is made to the section lift-curve 
slope. Thus, the value of m 0 to be used is the two-dimensional'value m Q 
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for sections normal to the elastic axis multiplied by a suitable correction 
factor. In reference 2, Shomick, in computing the divergence speed of 
unswept wings by strip theory, applies the aspect-ratio correction 


m e 



(V) 


For the case of an unswept uniform wing with 



this 


yields a divergence speed that differs by less than 1 percent 
result calculated by Hildebrand and Reissner (reference 3) by 
theory. (The parameter (3 calculated by lifting—line theory 


may be compared directly with ap 



of the present paper 


assumption 

from the 
lifting-line 
in reference 

.) 


Equation (4) may be extended to swept wings by means of the reasoning 
used in reference 6. If the air forces are considered to act in planes 
perpendicular to the elastic axis instead of parallel to the plane of 
symmetry, the following relation is obtained 


3 


* 




me 


= “o 


A 

A + 4 cos A 


(4a) 


This approx ima tion for finite—span effect may be used for subsonic and 
subcritical Mach numbers. For supersonic and supercritical Mach numbers, 
it is inapplicable; no span correction is available at present for these 
speed ranges. Neglect of the finite—span effects, however, will always 
tend to give conservative results. 


Location of the Effective Root 

If the results of the preceding analysis are applied to an actual 
wing, it is necessary to assume that the wing is clamped along a line 
perpendicular to the elastic axis. From the data and the analysis 
presented in reference 7 it appears that the amount of wing twist can 
be estimated closely by assuming the wing to be clamped at a line 
through the intersection of the elastic axis and the side of the fuselage, 
if the carry—through structure is fairly rigid. 

The bending deflections are estimated in reference 7 by considering 
the wing supported flexibly at a line through the innermost point of the 
trailing edge (in the case of a sweptback wing) or of the leading edge 
(in the case of a sweptforward wing). This effect may be taken into 
account in the preceding analysis by modifying the root boundary conditions. 
If the carry—through structure is fairly rigid the resulting divergence 
speed is the same, for all practical purposes, as the value obtained by 
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considering the wing mounted rigidly at a line through the intersection of 
the elastic axis and the side of the fuselage. This line•is consequently 
considered the effective root of an actual wing, as shown in figure l(c). 


APPLICATION OF THE THEORETICAL RESULTS 
Selection of the Aerodynamic Parameters 


Both the section lift-curve slope ihq and the aerodynamic center 
(and consequently the distance from the aerodynamic center to the elastic 
axis, ej) vary with Mach number and, in the transonic range at least, 
with the angle of attack. The parameters should be chosen at the angle 
of attack near the design value which yields the most conservative results. 
The choice of the Mach number depends on the purpose for which the 
divergence speed or dynamic pressure is calculated. 

If the divergence speed is calculated for its own sake, the parameters 
should be chosen at the Mach number corresponding to the divergence speed; 
resort to a trial and error procedure may therefore be necessary. If, on 
the other hand, it is desired to calculate the divergence speed or dynamic 
pressure as a reference value for one of the other aeroelastic phenomena, 
the aerodynamic parameters should be chosen at the Mach number of interest 
for the particular phenomenon. The dynamic pressure calculated in this 
manner will not usually be the true value but appears to be the divergence 
value only at the given Mach number and will vary with Mach number. If 
the variation of this reference divergence dynamic pressure qp is 

plotted against Mach number and the actual dynamic pressure q is plotted 
on the same graph, the intersections of the two lines will determine the 
true values of the divergence Mach number and dynamic pressure; this 
procedure therefore constitutes another way of calculating the divergence 
speed. Such a plot is shown in figure 5 for a straight and a sweptforward 
wing designed for high-speed flight. 

The values of the aerodynamic parameters are best obtained from 
experimental section data for the lift and pitching moment of the given 
section at the Mach number of interest. If such data are not available, 
the section lift-curve slope may be approximated by the G-lauert— 

Prandtl and Ackeret relations in the subsonic and supersonic regions, 
respectively, and by an arbitrary constant value which depends on the 
critical Mach number in the transonic region, as indicated in figure 6. 

The aerodynamic center is located near the quarter-chord point at subsonic 
speeds and at 40 to 45 percent of the chord at supersonic speeds. Its 
location at transonic speeds is not well established; it may move forward 
of the quarter—chord point, however, before receding to the supersonic 
position. 
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Use of the Curves 

The aerodynamic parameters are selected for a given Mach number and 
the section lift-curve slope is corrected for finite span effects in the, 
manner previously indicated. With these parameters and the given geometric 
parameters the ratio dp/ap may be determined from equation ( 3 ). The 

value of ap or dp is then obtained from figure 4 and the divergence 
dynamic pressure calculated from either of the two following relations: 


ID » a D - 

( 5 ) 

mgC-pL-^ cos A 1 r 


, Mr 

- dp cot A 

(5a) 

n^CpL - 3 cos^A 



The divergence speed is given by the relation 

= (6) 

If the value of qp as calculated in this manner is negative the 
wing cannot diverge, since a negative dynamic pressure does not correspond 
to any real speed. A negative value, however, is still useful as a 
reference value. 

Unless the wing under consideration is solid or has geometrically 
similar cross sections along the span, it must be kept in mind that the 
actual divergence speed or dynamic pressure may be below that calculated 
by this method for reasons cited subsequently in the "DISCUSSION." 


Use of Approximate Formulas 

The relation between the values ap and dp may be approximated 
closely by a straight line. The agreement between the exact curves and 
curves representing the linear relation is seen in figures 2 and 3 in the 
case of a uniform wing. From the linear approximation an expression for 
the divergence dynamic pressure may be obtained of the following form: 

* 

q = (GJ)r f _L_\ 

D m e c r l3 cos 2 A V e l c r/ 


K- 


1 - Kq 


(GJ)r / L \ 
(EI) r v e i c ry 


tan A 


(7) 
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where the constants Kp and K 2 are given in the following table: 


\ 

% 

Kg 

0.2 

2.81 

0.6i4 

..5 

2.74 

.497 

1,0 

2.47 

• 390 

1.5 

j 

2.22 

.326 

L__ 


Equation ( 7 ) may "be used with almost the same accuracy as the curves 
of figure 4. This equation should be particularly convenient for calculating 
the stiffnesses required for a given value of the dynamic pressure at 
divergence. 


EXPERIMENTAL RESULTS 


As a check on the theory several divergence tests on uniform swept 
wings were made in the Langley 4.5-foot flutter research tunnel. In the 
first series of tests a thin plate of 24S—T aluminum alloy, 5 inches 
by 30 inches by 0.126 inch, was held in the tunnel essentially as 
shown in figure l(a); the angle of sweep was varied by means of the 
rotating root fixture. The experimental divergence dynamic pressure was 
taken to be the highest value at which the wing would remain in an 
undeflected position when the root section was at zero angle of attack 
relative to the true air stream. From the experimental data obtained, 
shown in table I, the divergence dynamic pressure is plotted against the 
angle of sweep in figure 7(a). The variation of divergence dynamic 
pressure with sweep obtained from the theory of the present paper is 
also shewn in figure 7(a). No experimental data have been obtained for 
sweepback; the wing fluttered, rather than diverged, at zero sweep, and 
flutter would undoubtedly be critical at any angle of sweepback. However, 
the agreement between the present theory and experiment for angles of 
sweepforward up to 4-0° i 3 excellent. 

A similar series of tests was run for another rectangular aluminum 
wing, 4 inches by 24 Inches by O.O 977 inch. (See table I.) The results 
given in figure 7 (b) also show good agreement between theory and experiment. 

A third series of tests was run on flat—plate aluminum models, 

5 inches by 30 inches by 0.125 inch, that were not rectangular, but were 
clamped at the root, and cut at the tip parallel to the wind stream 
(fig. 1(c)). A separate model was used for each angle of sweep, but 
the length along the leading edge was held constant at 30 inches. These 
wing plan forms do not correspond to those assumed for the theory, but 
"they do represent more closely conventional swept^wing plan forms. The 
experimental results are presented in table I and plotted in figure 7(c). 

The theoretical variation of the divergence dynamic, pressure with sweep 





obtained, by using the effective root discussed previously (root A of fig. 7(c)) 
is shown by the solid—line curve. The variation obtained from theory if 
the effective root is considered to be at the line through the innermost 
point of the leading edge (root B of fig. 7 (c)) is shown by the dash-line 
curve. It appears that the use of root A yields conservative results for 
all sweep angles of practical concern, but that the results obtained by 
using root B are on the average in somewhat better agreement with the test 
data for low and medium angles of sweepforward. 

The curves of figure 7 are based on calculated values of (GJ) r 
and (El) r . The section aerodynamic center was assumed to be at 
the c/4 point, so that e^ is O. 25 . A section lift slope of 2it was 

assumed and corrected for span effects as shown in equation ( 4 a). 
Compressibility effects were neglected in all comparisons, since the 
test Mach numbers were too low to warrant a correction. 


DISCUSSION 


The aerodynamic and structural assumptions made impose limitations 
on the accuracy of the a n alysis. The aerodynamic assumptions are 
concerned with the magnitude of the deflections involved as well as the 
treatment of the aerodynamic induction and compressibility effects. The 
assumption of small deflections made in obtaining the effective velocity 
component V A (see appendix, equation (A5a)) yields results which are 

too low at high angles of sweep. The air forces are consequently 
underestimated and the divergence speed is overestimated. The correction 
for aerodynamic induction is only approximate; it is simple to make, 
however, and yields good experimental agreement. For most airplanes the 
divergence speed lies In the transonic or supersonic region, where the 
induction effects are greatly diminished, so that no correction need 
ordinarily be applied. The manner In which compressibility is taken into 
account in the analysis has not been checked e^erimentally and is 
consequently somewhat uncertain. 


The location of the aerodynamic center may not be known very 
accurately; fortunately, however, the divergence speed of a swept wing 
is not very sensitive to the aerodynamic-center location at large angles 
of sweep, as may be seen from figure 4(c). The change in d D , and hence 
the change in divergence speed, will be found to be fairly small even for 
rather large changes in the value of e;L produced by changes in the 

aerodynamic-center location or elastic-axis location. At low angles of 
sweep the effect of the location of the aerodynamic center relative to 
the elastic axis is much more pronounced; movement of the elastic axis 
forward or the aerodynamic center rearward will tend to raise the 
divergence speed. If the aerodynamic center is behind the elastic axis, 
negative values will usually be obtained for the dynamic pressure at 
divergence, so that divergence will be impossible. 


The most important of the structural assumptions involves the existence 
of a straight elastic axis, which permits tending and torsion to he treated 
independently• This assumption appears warranted in most cases. Although 
the exact location of this axis is often difficult to ascertain, only an 
approximate value is needed in view of the foregoing consideration, at 
least for large angles of sweep. Near the wing root, where the elastic— 
axis location is most -uncertain, -swept wings are so stiff that the location 
is immaterial. 

The assumption of the effective root at the location indicated in 
figure l(c) (root A of fig. 7(c)) leads to conservative results. This 
assumption appears, therefore, to he preferable to the assumption of the 
root at a point farther outboard, such as root B of figure 7(c). On the 
other hand, the assumption of the effective root B leads to results which, 
at least in the case of the flat—plate models tested, agree somewhat more 
closely with experiments for angles of sweepforward up to about 45°. 

The insensitivity of the divergence speed to the relative location 
of the aerodynamic center and the elastic axis indicates that the effects 
of bending predominate at large angles of sweep (actually at large values 
of dp/ap). Even at low values of sweep, bending is quite significant; 
the divergence phenomenon of swept wings should therefore be referred 
to as "bending—torsion" divergence rather than "torsional" divergence 
as in the case of unswept wings. 

The stiffness variation used in the analysis for tapered wings 
(El and GJ vary as c^) is realized for wings with geometrically 
similar cross sections; it is obtained for solid wings and closely 
approximates that of actual wings with a taper'ratio of the order of 0.2. 

For higher taper ratios, actual wings are more flexible at the tip than 
the fourth-power variation would dictate. In order to investigate this 
effect, the parameter ap has been plotted against taper ratio in figure 8 

for unswept wings with stiffness variations dictated by constant bending- 
stress levels. The assumption has been made that the torsional stiffness 
is proportional to the flexural stiffness which, in turn, is based on a 
load distribution given by strip theory. The computations for the 
divergence speed were performed by a method similar to that of reference 2. 
Figure 8 indicates that the divergence speed of wings with this type of 
stiffness distribution is in general lower than that obtained for wings 
the stiffness of which follows the quartic variation. The divergence speed 
of actual wings may be expected to lie between the two curves of figure 8 
and will in general be lower than that for wings the stiffness of which 
follows the fourth-power variation. This point must be considered when 
figure b is used; the results furnished by the curves of figure b may be 
somewhat unconservative for actual wings. 

The effect of sweep may be seen from figures 3* b, and 7- Figures 3 
and b show that as the sweep angle (and hence the parameter dp/ap) is 
increased negatively, that is, toward increased sweepforward, the 
parameter ap (and hence the divergence speed) decreases rapidly, all 


other physical parameters remaining the same. Beyond a certain angle of 
sweepforward the cos A term in the parameter a D becomes dominant 

and tends to increase the divergence speed again as seen in figure 7 . 
Sweptback wings with a value of dp/ap larger than about 2 and positive 
values of ep, which normally correspond to subsonic speeds, cannot 
diverge; nor can any sweptback or unswept wing diverge for negative values 
o ep which may exist at supersonic speeds. In either case a negative 
divergence dynamic pressure is obtained. This fact does not preclude 
the possibility of flutter in these cases, however. 

The effect of taper on the divergence of swept wings of the assumed 
stiffness variation (GJ and El varying as c^) is seen from figure 4. 

For positive values of e ± (subsonic case) and positive, or small negative 

values of the sweep angle, conventional taper increases the divergence 
speed; for all other configurations it decreases the divergence speed. 

The effects of inverse taper are the opposite of those of conventional 
taper. These considerations must be modified in the case of actual 
vings, because of their deviation from the assumed stiffness distribution 
so that either the curves of figure 4 must be used with some degree of ’ 
conservatism or more refined analyses must be resorted to. 

Although the low—speed wind-tunnel tests for divergence have been 
performed only on uniform wings, they serve to corroborate the theoretical 
analysis for uniform wings and hence, indirectly, that for tapered wings, 
since the assumptions are the same in either case. The rapid decrease 
of divergence speed.with increase in sweepforward agrees with the predicted 
variation both qualitatively and quantitatively (up to about 40° sweep- 
orward). This agreement indicates that the assumptions made concerning 
the structural and aerodynamic behavior (at low speeds) are justified. 

At values of sweepforward above about 40° the observed increase in 
divergence speed falls short of the predicted increase. This observation 
agrees qualitatively with the statement made previously that the analysis 
would tqnd to overestimate the divergence speed for large angles of sweep. 

It is difficult to make a direct comparison with the results of the 
simpiified approach of reference 5 , but it is quite apparent that since 
the method of this report takes account of the spanwise variation of 
bending and twisting it should furnish more reliable results. If the 
first natural bending and torsiop frequencies of a uniform cantilever 

"bco 

beam are substituted in the expression for —- of reference 5 a 

Tp * 

relation may be obtained for q D which will have both the same form as 
equation ( 7 ) and the same value of K ± if a proper finite-span correction 
is applied; the value of K^, however, will be 0.299 insteiad of O. 39 O, 

as determined by the theory of this paper for the uniform wing (X = l!o) 
so that the effect of sweep will be underestimated somewhat. This "' 

discrepancy indicates that, if the results of a semirigid analysis are 
applied to actual swept wings, a certain amount of caution must be exercised 


C ONCLIIDLNG- REMARKS 


On the basis of certain assumptions theoretical results have been 
obtained for the divergence of uniform and tapered wings. These results 
are presented in the form of charts and approximate formulas suitable 
for obtaining quick estimates of the divergence dynamic pressure or speed. 
A limited number of low-speed wind-tunnel tests on untapered models give 
good agreement with the theoretical results for angles of sweepforward 
up to k0° and appear to justify the assumptions made in the analysis. 

The results indicate the following conclusions: 

1. The divergence speed drops rapidly as sweepforward increases up 
to about ^0°. Wings with sweepback beyond a fairly low value cannot 
diverge. 


2 . Moving the elastic axis forward raises the divergence speed 
appreciably for low angles of sweep but has less effect at higher sweep 
angles. Wings with the elastic axis forward of the aerodynamic center 
(supersonic case) can only diverge for moderate or large angles of 
sweepforward. 

3. For most practical cases, the effect of conventional taper is to 
increase the divergence speed of essentially unswept wings and to decrease 
the divergence speed of wings with moderate and large angles of sweep if 
the stiffness varies as the chord to the fourth power. For stiffness 
variations closer to the ones which are obtained for actual wings, this 
effect may not be observed. In order to obtain accurate results in these 
cases a more refined analysis must be resorted to. 
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APEERDIX 


ANALYSIS OF THE D3YERGEHCE OF SWEPT WHJGS 


, Assumptions 

The analysis involves the following limitations and. assumptions! 

(a) Aerodynamic induction is taken into account only insofar as 

an over-all correction is applied to strip theory 

(b) Aerodynamic as well as elastic forces are based on the 

assumption of small deflections 

(c) The wing is clamped at the root perpendicular to a straight 

elastic axis (see fig. 1 ), and all deformations are 
considered to be given by the elementary theories of 
bending and of torsion about the elastic axis 


Aerodynamic Forces 

In keeping with assumptions (a) and (b), the force per unit width 
on a wing section in a plane perpendicular to the elastic axis resulting 
from bending and twisting deformations is given by 

i = *y* e fv A 2 c (ai) 

where <x e and are the effective angle of attack and the component 

of the free-etream velocity in the plane perpendicular to the e las tic 
axis, respectively. If the free—stream velocity is resolved into three 
components, one along the local tangent to the elastic axis, one parallel 
to the chord, and one perpendicular to the other two, the effective 
angle of attack may be obtained as the ratio of the third component to 
the second and the effective velocity as the vector sum of the two 
components. Thus, 


and 


tan a e 


sin cp cos A — cos cp sin f sin A 
cos cp cos A + sin cp sin f sin A 


V - vS jjcos cp C03 A + sin cp sin f sin A) 2 

+ (sin "j cos A — cos cp sin T sin/.) 2 ] 

The most convenient way of obtaining these relations is 
of vector analysis. 


(A2) 


(A3) 


probably by use 



(A4) 

(A5) 

(A5a) 


For small values of 9 and r these relations simplify to 

oc e = 9 — r tan A 
= v cos a 
= V cos A 

so xhat, for small deflections, equation (Al) can he expressed as follows 

2 = (<? - r tan A) 1 % ^V 2 cos3dc (A 6 ) 

Tlie moment of this force about the elastic axis is assumed to be 
given by the product Ze^c. The presence of a moment at zero lift due 

to camber will not affect the results in any manner, since the divergence 
speed depends only on the rates at which the aerodynamic and structural 
forces increase with the deformations. 


Differential Equations 

The differential equations of tending and torsion of the wing 
referred to the elastic axis are 


$ ( EI §) - 1 <^> 

m 

d7 ( GJ §) = “ Ze i c (A8) 


respectively. Substituting the small—deflection expression for Z 
(equation (A 6 )) into equations (A 7 ) and (a 8 ) yields the two simultaneous 
differential equations: 


~r (EI Uf = £v 2 cos 2 Am e c (9 — T tan A) (Ay) 

<ty 2 \ d 

h ( GJ df) = ~2 72 cos2Am e e l c2 & ~ F tan A) (A10) 


These equations are subject to the following boundary conditions: 


Zero local dihedral and twist at the root, 

r (0) = 0 


(All) 
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Zero moment, torque, and shear at the free end. 



(A12) 



0 



(A12a) 


(A12h) 


Equations (A9) and (A10) may he solved by numerical methods for 
any arbitrary stiffness and chord variation. In the case of untapered 
wings and linearly tapered wings with the stiffnesses varying as the 
chord to the fourth power, they can he solved directly, as shown in the 
following sections. 


Solution for Uniform Wings 

If the tending stiffness, the torsional stiffness, and the chord 
of the wing have constant values of (El) r , (GJ) r , an& c r> respectively, 
along the wing span, equations (A9) and (A10) become 

F 111 tan A = d(cp — F tan A) (A13) 

cp»* = —a(cp — T tan A) (Al4) 

where the differentiation denoted hy the prime is with respect to r\ = 
and the two dimensionless parameters a and d are defined hy 


a 


q cos^rr^eqCp L 


2 t 2 


(GJ) r 


(A15) 


q cos^VnigCpLStaxi A 

d = --- 

(El) r 


(A16) 


Differentiating equation (Al4) once and combining it with equation (A13) 
yields the single differential equation 

cxe* 11 + aa e * + da e = 0 . (A17) 


where a e - rp — F tan A. 


b-»H 


The general solution of this equation is 



i=l 

where the r^ f s are the roots of the characteristic equation 

r^ + ar + d = 0 


(A18) 

(A19) 


and. the Aj_ f s are arbitrary constants. 

The boundary conditions for equation (A17) are: 

From equations (All) and (Alla) 

Oe(0) = 0 (A20) 

From equations (A12) and. (A12a) 

a e *(l) = 0 (A20a) 

From equations (Al4) and (A12b) 

a e **(l)+aae(l) =0 (A20b) 

Substituting equation (Al8) into these boundary conditions gives 



The condition for divergence is that there be a nonvanishing 
solution for the Aj_*s; that is* a solution for which oc^ is not 

zero along the entire span. Therefore, one or more of the Aj/s must 
be different from zero (see equation (Al8)). Hence,'the determinant of 
the coefficients of the Aq f s in equations (A21) must vanish. Thus, 
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r l 2 + a ) e ' 


(r2^ + a)e r ^ 


= 0 (A22) 


To 1 - + a 


>)e r 3 


Critical combinations ap and d D of the parameters a and d 
are obtained if the combinations give rise to roots rand r^ 

for equation (A19) such that equation (A22) or its expanded equivalent 
is satisfied. The curve of ap against dp in f igure 2 constitutes a 
plot of these critical combinations; points on the curve vere computed 
by assuming values of ap and solving for dp by trial. 


Solution for Tapered Wings 

For tapered wings the chord varies linearly and the bending and 
torsion stiffnesses are assumed to vary as the fourth power of the chord 
Then, 


c = kc- 


where 


k = 1 - (1 - X)tj 

and X is the taper ratio c-^/c r . Furthermore^ 

GJ = (GJ) r k 4 

EE = (EI)^ 

The differential equations (A9) and (A10) then become 
tan A[k 3 r m + 8k 2 r u + 12kT f - d^rj + d^tp = 0 
jk 2 cn n + 4ktp' + arjrcp — arpF tan A = 0 


(A23) 


(A24) 


where 


(i - x ) 2 


a - 4 ) 3 


• • 
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The differentiations denoted by the primes in equations (A23) and (A2k) 
are vith respect to k rather than r). This procedure places the 
differential equations in the form of the Euler (or Cauchy) equations, 
which are easily tractable mathematically. 

Differentiating equation (A24) once, multiplying it by k, and 
combining the result with equations (A23) and (A24) yields a single 
differential equation in 


k 3 a e »*» + 8k 2 a e "+ (12 + a^ka-g* + (2aqi - dj,)^ = 0 


The solution of this equation is 


% = BikS1 
i =1 


(A25) 


(A2 6) 


vhere the Sj/s are the roots of 

s(s - l)(s - 2) + 8s(s - 1) + (12 + aqOs + (28^^ - dj) = 0 
or 

s 3 + 5^ + (6 + 8rji)s + (2ap — djj) = 0 

and the Bj/s are arbitrary constants. 

The boundary conditions are: 

From equations (All) and (Alla) 

a e (l) = 0 

From equations (A12) and (A12a) 

a e *(X) = 0 

From equations (A24), (A12a) and (A12b) 

* (X) + aQnCng(X) = 0 

Substituting equation (A26) into the boundary conditions gives 


B 1 = 0 


(A27) 


(A28) 

(A28a) 

(A28b) 


i=l 


Y~_ s i B i^- S1 = ° 

i=l 


3 


). _ | s i( s i - 1) + a T BjX bl = 0 

i=i 


(A29) 


I 
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Setting the determinant of the coefficients of the- B^'s in 
equation (A29) equal to zero yields, after some simplification. 


1 

1 

1 



X S 1 

s 2 X S2 

s 3 X s 3 . 

= 0 

(A30) 

• ar)*. 81 

^sg 2 + a T^ S ^ 

(s^ 2 + a r)l S 3 




+ 


Equation (A30) is very similar to the corresponding determinantal 
equation (A22) for untapered wings. For particular values of the taper 
ratio \, it determines critical combinations of a and d. Calcu¬ 
lations have been carried out for taper ratios of 0.2, 0.5, and 1.5, and 
the results are given by the curves of ap against dp/ap and dp 
against ap/dp in figures 4(a), 4(b), and 4(c), together with the 
results of the analysis for the uniform wing (taper ratio 1.0). 


m 
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TABLE I.- RESULTS OF DIVERGENCE TESTS 


-A 

(deg) 

<lD 

(experiment) 

(lb/sq ft) 

(calculated.) 

(lb/sq ft) 


Series 1 



5 in. by 30 in. "by 0.126 in.; square tips; e^_ = 0 . 25 ; 


(GJ) r = 13,330 lb-in. 2 ; (El) r = 8,830 lt^in . 2 

1 

_1 

0 

a 112.80 

178.6 

5.0 

83.20 

80.7 

14.7 

44.21 

40.9 

30.0 

26.50 

27.0 

45.0 

24.41 

26.2 

55-9 

24.84 

31.1 

55-9 

24.95 

31.1 

63.2 

26.65 

39-7 


Series 2 



4 in. "by 24 in. by 0.0977 in.; square tips; e-j_ = 0 . 25 ; 


(GJ) r = 4,980 lb-in. 2 ; (El) r = 3,300 lt^in . 2 


0 

a 106.90 

163.1 

5.0 

76.99 

73.8 

14.7 

43.10 

37.4 

30.0 

26.51 

24.7 

45.0 

22.44 

23.8 

60.0 

23.36 

32.0 

69.6 

23.74 

50.5 


Series 3 



5 in. by 30 in. by 0.125 in.; tips parallel to plane 


of symmetry; e^ = 0 . 25 ; 




(GJ) r = 13,050 lb-in. 2 ; (El) r = 8,560 lb-in. 




Root A 

Root B 

5.6 

79.87 

74.0 

74.8 

15.5 

41.80 

38.3 

40.8 

30.1 

31.72 

26.5 

30.3 

45.2 

28.90 

25.3 

32.8 

b 59.7 

c 38 .l 8 

33.7 

53.1 

74.6 

54.62 

79.3 

232.4 


a Flutter. v '-.nac a - 

^Model thickness, 0.123 inch. 
c Corrected for discrepency in thickness. 














(a) Assumed constant-chord (b) Assumed tapered wing, 

wing. 



Figure 1.- Definitions of geometrical parameters. 















Figure 3.- The relation between the parameters a and d at divergence for a uniform wing. 
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Figure 4.- Divergence of swept wings. 
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(b) Positive values of d D /a D . 
Figure 4. - C ontinued. 
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Mach number, M 

Figure 5.- Divergence of a straight and a sweptforward wing. 
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Taper ratio, \ 

Figure 8.- Effect of taper on the divergence of unswept wings. 





